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We compare the model power spectrum, computed based on perturbation theory with the power 
spectrum of luminous red galaxies (LRG) measured from the Sloan Digital Sky Survey Data Release 
7 catalog, assuming a fiat, cold dark matter-dominated cosmology. The model includes the effects of 
massive neutrinos, nonlinear matter clustering and nonlinear, scale-dependent galaxy bias in a self- 
consistent manner. We first test the accuracy of perturbation theory model by comparing the model 
predictions with the halo power spectrum in real- and redshift-space, measured from 70 simulation 
realizations for a cold dark matter model without massive neutrinos. We show that the perturbation 
theory model with bias parameters being properly adjusted can fairly well reproduce the simulation 
results. As a result the best-fit parameters obtained from the hypothetical parameter fitting recover, 
within statistical uncertainties, the input cosmological parameters in simulations, including an upper 
bound on neutrino mass, if the power spectrum information up to fc ~ 0.15 /iMpc~^ is used. 
However, for the redshift-space power spectrum, the best-fit cosmological parameters show a sizable 
bias from the input values if using the information up to fc ~ 0.2 /iMpc~^, probably due to nonlinear 
redshift distortion effect. Given these tests, we decided, as a conservative choice, to use the LRG 
power spectrum up to fc = 0.1 /iMpc~^ in order to minimize possible unknown nonlinearity effects. 
In combination with the recent results from Wilkinson Microwave Background Anisotropy Probe 
(WMAP), we derive a robust upper-bound on the sum of neutrino masses, given as < 0.81 eV 

(95% C.L.), marginalized over other parameters including nonlinear bias parameters and dark energy 
equation of state parameter. The upper bound is only slightly improved to ^ < 0.80 eV if 
including the LRG spectrum up to fc = 0.2 /iMpc~^, due to severe parameter degeneracies, though 
the constraint may be biased as discussed above. The neutrino mass limit is improved by a factor 
of 1.85 compared to the limit from the WMAP5 alone, X^m,y < 1.5 eV. 

PACS numbers: 98.80.Es,14.60.Pq,98.65.Dx 



I. INTRODUCTION 



Combining the cosmic microwave background (CMB) 
with large-scale structure probes provides a powerful 
means of constraining the sum of neutrino masses [1, 2]. 
Massive neutrinos imprint a characteristic suppression 
in the clustering of galaxies at scales below the free- 
streaming scale of neutrinos in a cold dark matter (CDM) 
dominated structure formation scenario. In particular, 
for neutrino masses of '^O.leV inferred from terrestrial 
experiments, a wide-field galaxy redshift survey can di- 
rectly probe the scales comparable with the neutrino 
free-streaming scale ^^lOOMpc, which is incidentally close 
to the baryonic acoustic oscillation scales. The existing 
galaxy surveys have provided a stringent upper limit on 
the neutrino mass ([3-5]; see Reid et al. 2010 [6] for 
the most recent study, hereafter RIO). However, all the 
previous studies employed a rather empirical approach 
to model the nonlinear effects in galaxy clustering such 
as the Qni-model [4] or the method to use polynomial 
functions of wavenumbcrs with additional nuisance pa- 
rameters [6]. We also note that these modelings have 



been tested using mock catalogs without neutrino effects 
being taken into account. 

In order to derive a robust, reliable constraint on neu- 
trino masses from the observed galaxy distribution, an 
accurate modeling of galaxy clustering is clearly needed 
properly taking into account the effects of nonlinear clus- 
tering, redshift distortion and nonlinear, scale-dependent 
galaxy bias. Simulation based approach may be the most 
powerful method, however, such a study for a mixed dark 
matter model (CDM plus massive neutrinos) is still in de- 
veloping stages, especially for small neutrino mass scales 
of interest ^ 1 eV [7-9] . An analytical approach is com- 
plementary, and allows us to study the effect of massive 
neutrinos as a function of different cosmological models. 
Recently wc have developed the new analytical method 
to compute the nonlinear galaxy power spectrum based 
on the perturbation theory (PT) approach [10, II], where 
the effects of nonlinear clustering and nonlinear galaxy 
bias are included in a self-consistent manner within the 
PT framework. The PT-based model (see also [12-14]) is 
a natural extension of the well-established linear theory, 
and the validity has been extensively studied by com- 
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paring with A'^body simulations in a CDM model (e.g. 
[15-18]). 

In this paper, wc present the first application of the 
PT model to the power spectrum of luminous red galax- 
ies (LRGs) measured from the Data Release 7 catalog 
of the Sloan Digital Sky Survey (SDSS) in RIO. We 
then derive a robust constraint on neutrino masses, com- 
bined with the WMAP 5-year (WMAP5) data [19], in- 
cluding marginalization over the uncertainties of galaxy 
bias parameters, residual shot noise contribution, and 
dark energy parameters. We mention the recent study 
[20], where the PT-based model is compared to the SDSS 
main galaxies (not LRGs). 

II. MODELING OF NONLINEAR GALAXY 
POWER SPECTRUM 

A. PT model 

In our previous papers [10, 11] we developed a PT- 
based method for computing the nonlinear galaxy power 
spectrum in a mixed dark matter model: 

Pg{k;z)^bl [pNL(fc;z) + fo2Pb2(fc;z) + 6^Pb22(fc;z)]+A^. 

(1) 

Here bi and 62 are the linear and nonlinear bias param- 
eters and N denotes the residual shot noise parameter, 
which are derived by renormalizing the galaxy bias pa- 
rameters based on PT prescription [21]. The expressions 
for power spectra P^^^, Ph2 and Pb22 are given in [11]. 
Note that Pb2 > and Pb22 < over a range of rele- 
vant scales. The nonlinear matter power spectrum P^ 
is given as 

+ 2/,b/,pL,(fc) + /2pL(fc)^ (2) 

where the subscripts "cb" and 'V denote "CDM plus 
baryon" and "massive neutrinos," respectively, and P^^ 

and P^^ describe the perturbative corrections to the 
power spectrum at next-to-leading order. The coefficient 
fi is the mass fraction of each species relative to the 
present-day energy density of total matter, fimo: fv = 
^Ml^mG = E"^<'/(^mO^^ X 94.1eV) and /^b = 1 - 
The nonlinear galaxy power spectrum at a given redshift 
z [Eq. (1)] can be computed once the linear-order power 
spectra of CDM, baryon and neutrino perturbations at 
the same redshift z are given for an assumed cosmological 
model and the bias parameters 61, 62 and N are specified, 
as extensively studied in [11]. 

B. Testing PT model with simulations 

To assess the validity of the PT model [Eq.(l)] in esti- 
mating model parameters, we implement a hypothetical 
experiment: By fitting the PT model to the halo power 



spectrum measured from TV-body simulations, we address 
whether the cosmological parameters assumed in TV-body 
simulations can be properly recovered. We used 70 simu- 
lation realizations each of which is carried out with 512'^ 
A^-body particles and volume of 1 h~^Gpc^, comparable 
with the volume covered by the SDSS survey (the A''-body 
simulations are kindly provided by Issha Kayo, and also 
see [22]). We created the halo catalogs from A^-body sim- 
ulation outputs at z = 0, based on the friend-of-friend 
method with the linking length of 6 = 0.2 (20% of the 
mean separation of A^-body particles). The catalog in 
each realization contains halos with masses greater than 
A'/niin — l0^^h~^MQ, where the mass threshold is de- 
termined such that the resulting number density of halos 
becomes fihaio — 3 x lO^^Ti'^Mpc"'^, comparable with that 
of the SDSS LRGs. 

Figure 1 shows the halo power spectra in real- and red- 
shift space, measured from the 70 simulation realizations 
above. The redshift-space power spectrum shown here 
is the monopole spectra, i.e. obtained by azimuthally 
averaging the 2D power spectrum in redshift space over 
circular annulus of a given radius of fc, where the line- 
of-sight direction is taken as the z-direction of each sim- 
ulation box. The filled circle at each k bin shows the 
mean band power computed from the 70 realizations, 
while the error bar is the scatter corresponding to the sta- 
tistical measurement uncertainty for a survey volume of 
1 [ft.~^Gpc]'^. The upper and lower panels show the real- 
and redshift-space power spectra, respectively, where the 
redshift-space power spectrum is affected by redshift dis- 
tortion effect due to peculiar velocities of halos. The 
redshift distortion effect considered here arises from the 
bulk motion of halos, because the peculiar velocity of 
each halo in simulation is defined by the mean of veloc- 
ities of member A'^body particles in the halo and there- 
fore the internal virial motions are averaged out. The 
linear theory predicts that the redshift distortion effect 
by the bulk motion causes only an overall shift in the 
power spectrum amplitude [23] , and does not change the 
shape of power spectrum. However, as can be found from 
the lower panel, the redshift-space halo spectrum, which 
is normalized so as to match the real-space power spec- 
trum amplitude at small fc, shows a scale-dependent en- 
hancement in the amplitudes compared to the real-space 
spectrum. This implies nonlinear redshift distortion due 
to the halo velocity field. (Also see [24] for discussion 
on velocity bias of halos.) Recently, the authors of [18] 
developed a model to compute the nonlinear, redshift- 
space power spectrum taking into account nonlinearity 
effects such as mass clustering and redshift distortion, 
based on perturbation theory. Interestingly, this work 
showed that such a scale-dependent enhancement in the 
redshift-space power spectrum amplitude may be caused 
by the nonlinear peculiar velocity, which qualitatively ex- 
plains the results shown in Fig. 1. However, a further 
exploration of this effect is beyond the scope of this pa- 
per, and will be studied elsewhere. 

In Fig. 1 we compare the halo spectra, measured from 
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the simulations, with the model predictions computed 
based on linear theory or the PT model. For the PT pre- 
dictions, we show the mass power spectrum (thin solid 
curve) assuming the cosmological model assumed in the 
simulations, and also show the best-fit halo power spec- 
trum (bold solid curve), computed from Eq. (2), where 
the best-fit model parameters including bias parameters 
are obtained by fitting the PT model to the simulation 
spectrum up to fc = 0.1 /iMpc~^ (also see below for de- 
tails). The figure clearly demonstrates that the simula- 
tion halo spectra cannot be explained by either the linear 
theory or the PT model for mass power spectrum over 
a range of wavenumbers of interest. Although the lin- 
ear theory may appear to give a good fit to the sim- 
ulation result up to fc ~ 0.12 hMpc'^, a closer look 
reveals that the linear theory over-estimates the band 
powers at scales around k ~ 0.07 /iMpc~^, where the PT 
model gives a better fit. On the other hand, interestingly, 
the PT-based halo spectrum with bias parameters being 
properly adjusted fairly well reproduces the simulation 
results. The PT model can give a reasonably good fit up 
to /c ~ 0.15 /iMpc~^, but begins to increasingly deviate 
from the simulation results at the larger k due to stronger 
nonlinearity effects. 

Next let us move on to details of the hypothetical 
parameter fitting; can the PT model for halo power 
spectrum recover the cosmological parameters assumed 
in the simulations? In the parameter fitting, we em- 
ploy 6 parameters: 17bo/f^mo(= 0.172), n^Qh{= 0.174), 
X] rrii, (= eV) , and the parameters of galaxy bias and 
shot noise, bi, &2 and N. (The values in parentheseis 
are the input values assumed in the iV-body simulations.) 
Other parameters are fixed to the input values of A'^body 
simulations. Note that adding the neutrino mass as a 
free parameter causes an asymmetric effect on the model 
spectra due to the sharp limit ^ rui, > 0: it causes only 
a scale-dependent suppression, not an increase, in the 
power spectrum amplitudes. We imposed the Gaussian 
prior a{flbo)/^ho = 0.05, and assumed the residual shot 
noise N to be smaller than the power spectrum ampli- 
tudes over the range of wavenumbers used. Finally, we 
used the power spectrum information up to the maxi- 
mum wavcnumber fcmax = 0.1 hMpc~^ , motivated by the 
fact that the PT model stays fairly accurate down to the 
wavenumber fc < 0.1 /iMpc^^ for mass power spectrum 
at z = in a CDM model as carefully studied in [16, 17]. 

Figure 2 shows the resulting posterior distribution of 
each model parameter. It is found that the input values 
of fibo/^mOj ^moh ^nd nii, are well recovered within 
the 68% C.L. statistical uncertainties as denoted by the 
shaded regions. A closer look reveals that the best-fit 
values slightly deviate from the input values. The ori- 
gin of the offsets can be explained by the dotted curves, 
which show the posterior distribution obtained by fixing 
the neutrino mass to zero (i.e., the input value for the 
CDM simulations) . Adding neutrino masses as a free pa- 
rameter causes such a bias in the best-fit values of ftho, 
f2mO^ and bi. This bias direction in flho, ^moh and bi is 
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FIG. 1: Upper panel: The filled circles at each k bin show 
the mean halo power spectrum measured from 70 simulation 
realizations at z = (see text for details), while the error bar 
shows the statistical measurement uncertainty at the k bin 
for a simulation volume of 1 hT^Gpc?, roughly comparable 
with the SDSS survey volume. For illustrative purpose the 
halo spectrum is divided by the no-wiggle, linear power spec- 
trum, multiplied by the linear halo bias squared, h\P^^^^{k) 
(6i = 1.66). For comparison, the thin-dotted and -solid curves 
show the linear-theory and PT predictions for "mass" power 
spectrum, respectively, for the cosmological model assumed 
in the simulations. The bold solid curve shows the best-fit 
PT model for halo power spectrum, computed from Eq. (2), 
where the best-fit model parameters including bias parame- 
ters are obtained by fitting the model predictions to the sim- 
ulation spectrum up to fe = 0.1 /iMpc~^ (see Fig. 2). Lower 
panel: Similar to the upper panel, but for redshift-space power 
spectrum {b\ — 1.81). The redshift-space power spectrum is 
modified by redshift distortion effect due to peculiar velocities 
of halos. For comparison, the circle points without error bars 
show the simulation halo spectrum in real space (the same as 
in the upper panel). 



found to all increase the power spectrum amplitudes so 
as to compensate a scale-dependent suppression caused 
by nonzero neutrino masses. Accordingly an upper limit 
on ^ is obtained due to the sharp cutoff ^ m^, > 0. 
It should also be noted that adding neutrino masses in- 
creases the marginalized error of bi, implying a strong 
degeneracy between bi and m^. 

A nonzero value of 62 is favored for the PT model to 
match the halo power spectrum or cquivalently a simple 
linear bias is disfavored even for fcmax = 0.1 hMpc^^, 
as also inferred from Fig. 1. The bimodal distribution 
of 62 is also apparent. For the favored values of 62, the 
term proportional to 62 in Eq. (1) is dominant over the 
term proportional to 62, and therefore both positive and 
negative values of 62 become acceptable. 

How is the parameter estimation changed if using the 
PT model up to higher fcmax than 0.1 ZiMpc"^, where 
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FIG. 2: Testing the perturbation theory (PT) based model 
with the halo power spectrum measured from TV-body simu- 
lations (70 realizations used). The solid curve in each panel 
is the posterior distribution of parameter, estimated by com- 
paring the PT model with the halo power spectrum up to 
the maximum wavenumber fcmax = 0.1 /iMpc~^. The input 
values of f2bo/f^mO and QmO, denoted by the vertical lines, 
are properly recovered within the statistical errors for the 
volume 1 h~^Gpc^, which is comparable with the SDSS vol- 
ume. For neutrino masses, which are not included in the 
A'^body simulations, an upper limit is derived. Nonzero val- 
ues of bias parameters (bi and 62) and shot noise parameter 
(A'^) are obtained, implying that the parameters are needed to 
describe the halo power spectrum. The dotted curves repre- 
sent the posterior distribution obtained by fixing the neutrino 
mass to zero, showing that the input values of the parameters 
flbo/^mo and fimo are correctly reproduced together with a 
tighter constraint on the linear bias parameter. 



the PT model ceases to be accurate, at least for the mass 
power spectrum [16, 17]? Some of the previous works 
sometimes attempted to use the power spectrum infor- 
mation up to such higher fc-range, motivated by the fact 
that the power spectrum of higher-fc modes contains a 
more constraining power of parameters. However, be- 
cause of complex nonlinearity effects, the best-fit param- 
eters derived from such high-fc information may be biased 
from the underlying true values. On the other hand, the 
PT model predicts a complex scale-dependent, nonlinear 
bias function function of cosmological model and 

bias parameters (see [11]). As implied in Fig. 1, the PT 
model can give a good fit to the simulation halo spec- 
trum at scales greater than k = 0.1 hMpc~-^, by adjusting 
the bias parameters. Therefore it is interesting to study 
whether or not fitting the PT model to the simulation 
halo spectrum up to the higher fc-values causes a bias in 
the best-fit parameters. 

Figure 3 shows the Icr marginalized errors on Sm^, 
^^bo/^mo and rimo^ as a function of the maximum 
wavenumber fcmax employed in the parameter fitting. 



First, let us focus on the results for the real-space power 
spectrum. Even though the PT model breaks down 
a.t k Si 0.1/iMpc~^ and over-estimates the "mass" power 
spectrum amplitudes at such high-fc range (see the up- 
per panel of Fig. 1), the best-fit parameters are found to 
recover the input values within the l-cr statistical uncer- 
tainties. It is also clear that the statistical errors of the 
parameter and the upper bound on Em^ arc improved 
at fcmax > 0.1 h^^Mpc~^ compared to our fiducial choice 
of fcmax = 0.1 /iMpc~^, due to a gain in the constraining 
power contained in the high fc-range. This may rcfiect 
that the PT model has more degrees of freedom to de- 
scribe the nonlinear halo power spectrum by adjusting 
the bias parameters, which may allow one to overcome 
the limitation of PT model for mass power spectrum. 

However, this is not the case for the redshift-space halo 
power spectrum. Again note that, if the halo power spec- 
trum measured from simulations is affected only by the 
Kaiser effect, the redshift distortion effect causes only 
an overall shift in the power spectrum amplitude, inde- 
pendent of fc, which can be absorbed by changing the 
linear bias parameter in the PT modeling. The fig- 
ure shows that the input parameters are recovered up 
to fcmax — 0.15 /iMpc~^, but the best- fit parameters at 
fcmax = 0.2 hMpc~^ show a sizable deviation, more than 
the statistical uncertainties, compared to the input val- 
ues for fimo^^ and 17bo/i^mO- This deviation implies that 
the residual nonlinear redshift distortion cannot be de- 
scribed by the PT model, even if changing the model 
parameters. 

Note that, on the contrary, at low fcmax ~0.05 hMpc~^, 
there are less statistical powers, giving larger uncertain- 
ties in model parameters. In addition, severe parameter 
degeneracies give only a very weak upper bound on neu- 
trino mass, and cause a bias in flmo^ as discussed above. 

Given the results in Figs. 1 and 2, we will use, as a 
conservative choice, the SDSS LRG power spectrum up 
to fc = 0.1 /iMpc~^ to compare with the PT model. For 
the LRG power spectrum measured in RIG, the nonlinear 
redshift distortion, known as the Fingers-of-God effect, is 
suppressed by clipping possible satellite LRGs. However, 
as we have shown, there may be a residual contamina- 
tion from the nonlinear redshift distortion effect. There- 
fore, in order to derive a robust constraint on neutrino 
mass, we will adopt fcmax = 0.1 /iMpc^^ for the following 
results, although we will also discuss how the neutrino 
mass constraint is changed by including the information 
beyond fc = 0.1 /iMpc^^, to be more comprehensive. 



III. RESULTS 

We now apply the PT model to the power spectrum 
of SDSS LRG samples in order to constrain neutrino 
mass. We use the halo power spectrum measured by 
RIO, where 104,337 halos were first reconstructed from 
the observed 110576 LRGs' distribution, and the angle- 
averaged redshift-space power spectrum was estimated 
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FIG. 3: The best-fit parameters and the marginahzed errors 
obtained by fitting the PT model with the halo spectrum 
up to a given maximum wavenumber fcmax, denoted in the 
horizontal axis. For each fcmax the left-side point with error 
bar shows the results for the real-space halo spectrum, while 
the right-side point shows the results for redshift-space halo 
spectrum. The horizontal dashed line denotes the input value 
of each parameter. 



based on the method in [25]. The halo power spectrum 
is less affected by the Fingcrs-of-God effect, because the 
contribution from satellite galaxies was eliminated in the 
measurement. Thus the halo power spectrum in RIO is 
appropriate to compare with the PT model. 

In estimating parameters, we combine the LRG power 
spectrum with the WMAP5 data. Note that our re- 
sults would remain almost unchanged even with the lat- 
est WMAP7 result [26] . We assume that the likelihood 
function of the LRG power spectrum is given as 



21n/:sDsscx 



(3) 



kij<k„ 



where is the difference between the measured and 
model power spectra at the i-th wavenumber bin ki, 
Ai = Phaio(fci) - P^^^{aki\p), with p being a set of 
model parameters (see below). Note that the effect of 
survey window function is properly taken into account 
in computing the model power spectrum following RIO. 
The matrix C is the covariance matrix for which we use 
the matrix provided in RIO, and is its inverse ma- 
trix. Note that we employ fc^ax = 0.1 hMpc~^ and as- 
sume the single redshift slice z = 0.35 for simplicity. The 
stretch factor "a" in the argument of the model power 
spectrum describes the cosmological distortion [27, 28]. 
This factor is given as a = D'^^{z)/Dv{z;p), where 
Dv{z) = [(1 + zfDA{zfz/H{z)Y^^ and D'^^ is the dis- 
tance for the reference cosmological model used in the 
LRG spectrum measurement. The likelihood for the 
joint analysis of WMAP5 plus SDSS is simply given as 
ln£ = ln£sDss + ln£wMAP- 

We include a fairly broad range of parameters that 
can cover variants of CDM cosmology such as models 
including massive neutrinos and dark energy equation of 
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FIG. 4: The parameter constraints obtained by comparing 
the PT model with the SDSS LRG power spectrum up to 
fcmax = 0.1 hMpc~^, in combination with the WMAPS con- 
straint, where we include 12 parameters given by Eq. (4). 
The upper panel shows the posterior distribution of neutrino 
masses, yielding the upper limit "^m,^ < 0.81 eV (95% C.L.), 
a factor 1.85 improvement over the limit J^m^ < 1.5 eV 
from the WMAP5 alone. The lower two panels show how 
the neutrino mass is degenerate with the dark energy equa- 
tion of state parameter w and the nonlinear bias parameter 
62, respectively, with 68% C.L. (dark shaded) and 95% C.L. 
(light shaded) regions. A nonzero 62 or equivalently a scale- 
dependent bias is favored at 68% C.L. Our results are com- 
pared with the results derived using the halo-model based 
method in RIO for the same maximum wavenumber cutoff 
fcmax = 0.1 hMpc-K 



state parameter. We vary 12 model parameters in total: 

p = {f„,Q.hoh^,^-DMoh'^,(^*,w,T,A'^,ns,Asz,bi,b2,N), 

(4) 

where flBMoh^ is the sum of CDM and massive neutri- 
nos: 9^ is the parameter to characterize the angular scale 
of CMB acoustic oscillations: r is the optical depth to 
the last scattering surface: and A^ are the parame- 
ters to specify the primordial power spectrum following 
the convention in [19]: Agz is the parameter to control 
a contamination of the Sunyaev-Zel'dovich effect to the 
CMB spectrum: w is the dark energy equation of state 
parameter. Note that the parameters r and Asz affect 
only the CMB information. We used the COSMOMC 
code [29] to explore parameter estimations in the multi- 
dimensional parameter space. 

The upper panel of Fig. 4 shows the marginalized 
error on neutrino masses. We obtain the upper limit 
Y^m^ < 0.81 eV (95% C.L.) for the SDSS Data Re- 
lease 7 plus WMAP5. This is a factor of 1.85 improve- 
ment compared to the limit derived from the WMAP5 
alone, 'Yrrii, < 1.5 eV. Our neutrino mass limit can 
be compared with the result derived using the method 
in RIO, where the empirical treatment based on halo 
model prescription was used to account for nonlinear. 
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FIG. 5: Comparing the best-fit PT model with the SDSS LRG 
spectrum, where the best-fit model is obtained fi:om the fitting 
to fcniax = 0.1 /iMpc~^. For illustrative clarity the power 
spectra are divided by the linear matter power spectrum for 
the best-fit cosmological model. For comparison, we also show 
the PT model, where the neutrino mass is changed to niv — 
0.81 eV, corresponding to the 95% C.L. upper bound in Fig. 4, 
but other parameters are kept fixed to the best-fit values. 



scale-dependent galaxy bias. Note that RIO employed 
kmnK = 0.2 /iMpc~^ and then derived an upper bound on 
the neutrino mass given as ^ < 0.63 eV. 

Fig. 5 shows that the best-fit PT model matches the 
measured LRG power spectrum well over a range of the 
working wavenumbers, k < 0.1 /iMpc~^. If the neu- 
trino mass is changed to the 95% C.L. upper bound, 
^ = 0.81 eV, but other parameters are kept fixed 
to the best-fit values, the model spectrum significantly 
underestimates the measured spectrum amplitudes at the 
small scales. Also note that the best-fit model rather con- 
tinues to match the measured spectrum beyond fcmax = 
0.1 hMpc^^. In fact, even if including the information 
up to /cmax = 0.2 /iMpc~^, the neutrino mass limit is 
only slightly changed to J2 "^i' — 0-^ reflecting less 
cosmological information at the higher wavenumbers due 
to severe degeneracies of cosmological parameters with 
nonlinear bias parameter and/or shot noise parameter. 

The lower two panels of Fig. 4 show how the neutrino 
mass constraint is degenerate with w and the nonlinear 
bias parameter &2- The marginalized constraint on w is 
-1.08 < u; < -0.79 (68% C.L.). While a change of w 
leads to a suppression in the power spectrum amplitudes. 
Fig. 4 shows that degeneracy between w and the neutrino 
mass is rather weaker than expected. This implies that 
the constraint on w comes mainly from the baryon acous- 
tic oscillation information as studied in [28] . Figure 4 also 
shows that a simple linear bias model with 62 = is disfa- 
vored at 68% C.L. That is, the nonlinear scale-dependent 
bias is needed to match the measured power spectrum, 
as can be found from Fig. 5. Similar to Fig. 2, bimodal 



structure of the constraint on 62 is found, implying that 
the term proportional to &| in Eq. (1) is dominant over 
other terms in the nonlinear power spectrum. 



IV. SUMMARY 

In this paper, we explore the robustness of the PT- 
based model to interpret the measured galaxy power 
spectrum, focusing on constraining the neutrino mass. 
The model successfully include the effects of nonlinear 
clustering and nonlinear, scale-dependent galaxy bias in 
a self-consistent manner within the PT framework. We 
have tested the accuracy of the PT model by compar- 
ing the model predictions with the halo power spectrum 
measured in the A-body simulation without massive neu- 
trinos. A careful and detailed comparison shows that 
the PT model can reproduce the simulated halo power 
spectrum and recover the cosmological parameters input 
in the simulations within statistical uncertainties, if the 
power spectrum is used up to fc ~ 0.15ft.Mpc~^. How- 
ever, in the case of the redshift-space power spectrum, 
the best-flt cosmological parameters show a biased esti- 
mation from the input values if the information up to 
k ~ 0.2/iMpc~^ is used. Thus, it is unclear to choose 
the maximum range of wavenumber but we decided to 
conservatively use the observed power spectrum up to 
fc = 0.1 /iMpc"^ in order to minimize possible unknown 
nonlinear systematic effects. 

Based on the test against the mock power spectra, we 
have applied this PT model to the SDSS LRG samples, 
and derived the neutrino mass limit ^ rrii, < 0.81 eV 
(95% C.L.). The parameter constraints including neu- 
trino masses would be further improved by including the 
redshift distortion measurement and/or the higher-order 
clustering information, which help to break the degenera- 
cies with galaxy bias parameters. On a theory side, the 
PT-based model needs to be further refined by including 
higher-order loop corrections and/or by calibrating the 
model with a suit of high-resolution simulations (see [18] 
for such a study) , although a careful treatment of massive 
neutrinos, nonlinear galaxy bias and redshift distortion 
is definitely needed. Once such refined models are avail- 
able, a more stringent constraint on neutrino masses can 
be obtained from high-precision measurements of galaxy 
clustering via on-going and future galaxy redshift sur- 
veys. We hope that this paper gives the first attempt to 
step in this direction. 
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